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ABSTRACT 

We have observed the most distant (= 0.829) cluster of galaxies in the Einstein 
Extended Medium Sensitivity Survey (EMSS), with the ASCA and ROSAT satellites 
We find an X-ray temperature of 12.31^ keV for this cluster, and the ROSAT map 
reveals significant substructure. The high temperature of MS1054-0321 is consistent with 
both its approximate velocity dispersion, based on the redshifts of 12 cluster members 
we have obtained at the Keck and the Canada- France- Hawaii telescopes, and with its 
weak lensing signature. The X-ray temperature of this cluster implies a virial mass 
^ 7 4 x 10 14 fi” 1 A/q, if the mean matter density in the universe equals the critical 
value (flo = 1), or larger if fi 0 < 1. Finding such a hot, massive cluster in the EMSS 
is extremely improbable if clusters grew from Gaussian perturbations in an fi 0 = 1 
universe. Combining the assumptions that D 0 = 1 and that the intial perturbations were 
Gaussian with the observed X-ray temperature function at low redshift, we show that this 
probability of this cluster occurring in the volume sampled by the EMSS is less than a few 
times 10~ 5 . Nor is MS1054-0321 the only hot cluster at high redshift; the only two other 
z > 0.5 EMSS clusters already observed with ASCA also have temperatures exceeding 
8 keV. Assuming again that the initial perturbations were Gaussian and fl 0 = 1, we find 
that each one is improbable at the < 10 -2 level. These observations, along with the 
fact that these luminosities and temperatures of the high-z clusters all agree with the 
low -2 Lx - Tx relation, argue strongly that fi 0 < 1. Otherwise, the initial perturbations 
must be non-Gaussian, if these clusters’ temperatures do indeed reflect their gravitational 
potentials. 


1. INTRODUCTION 

Clusters of galaxies occupy a special position 
among celestial self-gravitating structures, as they 
are the largest virialized objects in the universe. 
If all the cosmic structures we currently see grew 
hierarchically from small seed perturbations in the 
early universe, as many lines of evidence now suggest 
(e.g., Peebles 1993), then the largest clusters of 
galaxies are also the most recent virialized structures 
to have appeared on the cosmological scene. The 
most massive of all clusters, which represent the most 
extreme high-density peaks on scales ~ 10 15 A/q in 
the initial field of density perturbations, correspond 
to several-<7 events if the perturbations axe Gaussian. 

Because massive clusters represent the high- 
er tail of the distribution, their number density 
should rise dramatically as the overall amplitude 
of the perturbations grows. This characteristic of 
cluster formation renders massive cluster evolution 
particularly sensitive to the mean matter density of 
the universe, currently equal to fio* in units of the 
critical density 3// 2 /8 ttG. If Qq % 1. perturbations 
should still grow like (1 + z)” 1 , and evolution in 
the number density of distant clusters with redshift 
should be rapid. If fio < 1, perturbation growth 
has virtually stopped, and massive cluster evolution 
should be much more mild. 

Evolution in the distribution of cluster masses 
sensitively indicates fid but measuring the virial 
masses of clusters (M v j r ) is difficult, even at 


low redshifts. A much more straightforward 
observational task is to measure how the X-ray 
temperatures (Tx) of clusters depend on redshift 
and then to use a model to relate Tx to M vir (e.g, 
Oukbir & Blanchard 1992). The Advanced Satellite 
for Cosmology and Astrophysics (ASCA) (Tanaka, 
Inoue, & Holt 1994) now makes such tests possible 
because it can measure the temperatures of the most 
luminous X-ray clusters at high redshifts ( Yamashita 
1994; Furuzawa et al. 1994; Donahue 1996). Henry 
(1997) has recently published a temperature function 
for 10 massive clusters at z ~ 0.3, chosen from 
the Einstein Extended Medium Sensitivity Survey 
(EMSS) (Gioia et al. 1990; Henry et al. 1992). His 
comparision of the z « 0.3— 0.4 temperature function 
with the z « 0 temperature function suggests that 
fio ^ 0.5 and rules out fto = 1 with 99% confidence. 

In this paper we attempt to extract the maximum 
leverage out of this test for Q.q by measuring 
the temperature of the highest-redshift cluster 
in the EMSS, MS1054-0321 at z = 0.83. We 
present both an ASCA spectrum and a ROSAT 
HRI image of this very luminous cluster ( L x = 
5.6 x 10 44 h~ 2 ergs" 1 2-10 keV rest frame, L x = 
1.1 x 10 45 /i“ 2 ergs" 1 bolometric). For qo = 0.1, 
L x = 7.7 x 10 44 h“ 2 ergs'" 1 2-10 keV rest frame. 
Luppino & Kaiser (1997) provide a true color 
optical image and a weak lensing map. After 
presenting the observations in § 2, we examine 
how they constrain fio in § 3, finding that they 
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Figure 1 — We plot the X-ray spectra from each 
of the four ASCA telescopes for the cluster 
MS1054-0321. Each spectrum has over 1000 X-ray 
counts. (See Table 1.) The best-fit model is drawn 
with a solid line. The spectra we display were binned 
such that each energy bin contained a minimum of 
16 (SIS) or 25 (GIS) counts. 

intrinsic absorption component, to a 90% upper 
limit of Nh < 1-4 x 10 21 cm -2 . Thus, we fixed soft 
X-ray absorption at an HI value of 3.7 x 10 20 cm -2 
(Gioia et al. 1990). Each SIS spectrum was fit with 
its own normalization; the normalizations of the 
GIS spectra were constrained to be the same. The 
parameters varied to provide a fit to 4 spectra were 
therefore the 3 independent normalizations (1 for 
each SIS spectrum, 1 for both GIS spectra), a single 
emission-weighted temperature, and the metallicity 
(iron abundance). Examples of the binned spectra 
and best-fit model convolved with each telescope’s 
response matrix are found in Figure 1. The 
best-fit redshifted temperature to the binned-by- 
four spectrum was 12.4±2l3 keV with an upper 
limit on iron abundance of 0.22 solar, keeping 
the absorption fixed. If we used the SIS spectra 
regrouped to a minimum of 1G counts per bin, 
the best-fit temperature was virtually the same: 
12.3±2 1 keV with an upper limit on the iron 
abundance of < 0.25 solar. All uncertainties quoted 
are the 90% confidence levels for a 1-dimensional 
fit (Ax 2 = 2.70), and all fits have an acceptable 
reduced \ 2 of - 0.8 - 0.9. Figure 2 shows the two- 
parameter x 2 contours for the cluster metallicity and 
Tx in keV. 





kT(keV) 

Figure 2 — We plot the two-dimensional x 2 con ~ 
tours at 68.3%, 90% and 99% confidence levels 
(Ax 2 = 2.30, 4.61 and 9.21) for the cluster iron 
abundance in units of the solar abundance and 
temperature in keV. These contours represent the 
confidence contours of simultaneous fits to the SIS 
spectra and to GIS spectra binned such that each 
bin contained a minimum of 16 and 25 counts 
respectively. 

2.2. Optical Convboration 

The available optical data on this cluster cor- 
roborate the high temperature measured by ASCA. 
Weak lensing observations of MS 1054-0321 by Lup- 
pino & Kaiser (1997) show that its mass within 
0.5 /i*" 1 Mpc is (5- 11) x 10 14 /i“ l Mpc, depending 
on the redshift distribution of the lensed galaxies. 
In an isothermal potential, its mass distribution 
corresponds to an isotropic one-dimensional velocity 
dispersion ~ 1100 - 2200 km s' 1 , consistent with 
the value expected from the ASCA temperature: 
(kTx/^Tn p ) 1 ^ 2 = 1400±170 km s'* 1 (90% confidence 
limits). 

In addition, we have obtained the redshifts of 
twelve cluster members with the CFHT Multi- 
Object Spectrograph (MOS) on April 13, 1994 
and with the LRIS (Oke et al. 1995) at Keck on 
January 4, 1995. These are listed in Table 2. Using 
Timothy Beers 1 program ROSTAT, we calculated 
their mean redshift and velocity dispersion, finding 
z = 0.829 ± 0.008 and a velocity dispersion of 
1360 ± 450 km s” 1 (90% confidence interval). (The 
observed redshift dispersion has been corrected to 
first order for cosmological redshift by dividing by 
1 + z). Several alternative methods (Beers, Flynn, 
& Gebhardt 1990) give the same answer within the 
uncertainties. While the small numbers of redshifts 
do not constrain the velocity dispersion very well, 
the value we find is consistent with both the lensing 
results and the X-ray temperature. 
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Plate 1 — The central 3.75’ x3. 75’ of the cluster MS1054-0321. The optical image is a 14,400 second I-band 
image taken with the University of Hawaii 88-inch (Gioia & Luppino 1994). The ROSAT HRI image, a 
sum of all available HRI data for a total livetime of 121,590 seconds, was rebinned into 8" x 8" pixels 
and adaptively smoothed with three Gaussians with sigmas of 16, 14.3, and 12.6 arcseconds corresponding 
to counts in a pixel of < 10, 10-13.5 and > 13.5. The sigmas were chosen so that the products of the 
number of counts and a 2 were roughly constant. The background was 9.4164 counts per pixel. The net 
contour levels are 1, 2.8, 4.6, 6.4, 8.2, and 10 counts per pixel. These net count rates correspond to X-ray 
surface brightnesses of 1.4, 3.9, 11.6, and 14.0 x 10 -14 erg cm -2 s -1 arcmin -2 assuming that 1 HRI count 
~ 2.8 x 10 -n erg cm -2 . The galaxies are marked with identifications corresponding to Table 2. The central 
galaxy is at RA(2000)=10 56 59.9 and Dec(2000) = -03 37 37.3, as measured from our HST image of this 
cluster (Donahue et al., in preparation). 
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The asymmetric appearance of the cluster in 
these high resolution images warns us that the 
relationship between its mass and its temperature 
might not be as simple as a spherical, hydrostatic, 
isothermal model would predict. However, recent 
hydrodynamic simulations of cluster formation be- 
ginning with cosmological initial conditions indicate 
that the systematic errors on cluster masses derived 
using a scaling relation between mass and tempera- 
ture are relatively small, typically < 20% even when 
some substructure is present (Evrard, Metzler, & 
Navarro 1996). Higher resolution simulations of 
individual cluster mergers, whose initial conditions 
are not dictated by hierarchical structure formation, 
indicate that errors in the estimated mass can be 
somewhat larger (,$, 50%) after a recent merger, but 
these errors tend to yield underestimates of the clus- 
ter’s mass (Roettiger, Burns & Loken 1996; Schindler 
1996). The following discussion will therefore focus 
on the implications of the cluster temperature alone, 
assuming that the emission-weighted temperature is 
representative of the virial temperature. Mass errors 
at the 20% level do not affect the qualitative results 
of the ensuing discussion. 

3. CONSTRAINTS ON TEMPERATURE 
EVOLUTION 

This cluster, MS1054-0321, is one of the hottest 
clusters known. If fl 0 = 1» structure formation 
simulations show that we can estimate its virial 
mass by assuming that the mean density within 
its virialized region is ~ 200 times the critical 
density at the cluster’s redshift and that the 
cluster itself is isothermal (Evrard et al. 1996; 
Eke et al. 1996). The virial mass of MS1054-0321 
would then be approximately Af vir ~ [2 kTx/prn p 
(1 + z)] 3/2 (10/f o G) _1 ~ 7.4 x 10 14 h -1 Af© within 
r 2 oo = 1.5h~ l Mpc. (We note that this relation is 
equivalent to assuming (3 = 1 in a standard X-ray 
mass relation for an isothermal gas.) Such a massive 
cluster at such a high redshift severely challenges 
models of hierarchical structure formation with 

= 1 (Evrard 1989; Peebles, Daly, k Juszkiewicz 
1989; Donahue 1993). This section elucidates how 
severe those challenges are. 

Because the virial masses of clusters are much 
harder to measure than their temperatures, we 
will restrict this analysis to focus on temperature 
evolution, relying on T\ to be a surrogate for 
A/vir- The analysis will therefore depend on 
the temperature-mass relation we assume. In an 
Do = 1 universe, this relation is relatively simple: 
T x °c + z). If < It this relationship 

breaks down below z ~ fig 1 “ 1 3,8 the development 
of structure stagnates. 


-V. 1. Rarity of Hot Clusters at High-z 

In order to determine the rarity of a cluster as hot 
as MS 1054-0321 in an flo = 1 universe seeded with 
Gaussian perturbations, we will assume that the 
Press-Schechter formula (Press k Schechter 1974) 
adequately describes the mass function of virialized 
objects, an assumption borne out by numerical 
simulations (e.g., Lacey k Cole 1994, Eke et al. 
1996). In integral form, the comoving mass density 
of virialized objects with virial masses greater than 
Af is then 

p(> M) = poei(c(u c /V2) = ^ / e -1 * dx , 

V 77 K/s/2 

( 1 ) 

where po is the current mean matter density and 
v c is the critical threshold at which perturbations 
virialize. This threshold is v c = 6 c (t)/rr(M)D(t) t 
where &(M) describes the rms amplitude of linear 
perturbations on mass scale A/, D(t) describes their 
linear growth rate, and 6 c (t) is the linear overdensity 
at which perturbations virialize. If fio = 1, then 
v c oc (1 4- z) at a fixed mass scale Af . 

This integral form of the Press-Schechter formula 
makes its theoretical underpinnings somewhat more 
apparent than the more familiar differential form 
does. The formalism assumes that the original 
linear density perturbations from which structure 
grew, smoothed over a length scale corresponding 
to mass A/, follow a gaussian distribution with 
a dispersion proportional to rx(Af). As the 
perturbations grow, the most prominent of them 
collapse and virialize when their amplitudes exceed 
some threshold. The parameter v c expresses how 
many standard deviations away from the mean 
amplitude a positive density perturbation must lie 
to have collapsed and virialized by time t. The 
total mass density in virialized perturbations of 
mass Af or greater is then proportional to the 
integral over the tail of this Gaussian distribution 
from v c to infinity, which yields the complementary 
error function in the above expression if cr(A/) — ► 0 
as Af oo. Because the typical amplitudes of 
perturbations decline with increasing Af , the number 
density of virialized objects exceeding lhass Af is 
approximately Af). At several standard 

deviations away from the mean, virialized objects 
with masses significantly exceeding Af are much 
rarer than objects of mass Af , so this approximation 
for the mean number density of virialized clusters is 
good to a few tens of percent. 

We will now proceed to determine the expected 
number density of clusters like MS1054-0321 at high 
redshift as follows: 
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virialized mass density in such clusters is similar 
to the case of MS1054-0321, implying u c > 3.4 at 
z = 0 and v c > 5.2 at z — 0.54. For Qo = 1 
and Gaussian perturbations, the expected number 
density of > 8 keV clusters at z — 0.54 is then 
< 8.6 x 10" 11 h? Mpc -3 . Each of the two hot EMSS 
clusters at z — 0.54 would therefore be rarer than 
one in a hundred, with a joint improbability of order 
a few times 10” 4 , somewhat more probable than the 
chance occurrence of MS1054-0321. 

In summary, given the three hot high redshift 
clusters in the EMSS, the implied number density 
of clusters with Tx > 8 keV at z — 0.5 — 0.9 is 
~ 1.6 x 10” 8 /i 3 Mpc -3 . The probability of all three 
clusters appearing in the volume sampled by the 
EMSS in an Q = 1 Universe is ^ 10“ 6 if the initial 
density perturbations have a Gaussian distribution. 

Universes with Qo % 0.3 have no problem 
accomodating such hot clusters at these redshifts. 
Eke et al. (1996; see also Viana & Liddle 1996) show 
that the cluster temperature function evolves very 
little in an Qo = 0.3 universe with no cosmological 
constant and only modestly in a flat universe with 
the same Qo. The number density of > 8 keV 
clusters from Henry (1997) is ~ 5 x 10” 8 h 3 Mpc -3 , 
quite close to the value we find at z ^ 0.5. The 
statistics of hot clusters at z > 0.5 therefore add 
substantial weight to the growing body of evidence 
emerging from cluster evolution that Qo 22 0.3 — 0.5 
(e.g., Carlberg et al. 1997; Henry 1997; Bahcall, Fan, 
& Cen 1997). 

3.3. High-z L\~Tx Relation 

Cluster evolution from z 0.5 to the present 
should be modest at best in a universe with low 
Q 0 . Here we show that the relationship between 
the X-ray luminosity (Lx) and the X-ray tempera- 
ture of MS1054-0321 is similar to the relationship 
derived from lower- redshift clusters. Edge & Stewart 
(1991b) found that L x - (8.8 x 10 41 /i -2 ergs -1 ) 
(T x /I keV) 2 - 79±0 05 in the 2-10 keV band for low- 
redshift clusters, while Mushotzky & Scharf (1997) 
demonstrate that clusters observed at z > 0.14 
by ASCA, including MS0451-03 (Donahue 1996) 
and MS0016+16 (Yamashita 1994; Furuzawa et al. 
1994), do not deviate appreciably from this relation. 
Our observations indicate that Lx (2 — lOkeV) ^ 
7.7 xlO 44 h“ 2 ergs" 1 for MS1054-0321 in its own rest 
frame (qo = 0.1), implying an expected temperature 
of 11.3 keV, on the lower edge of our one-dimensional 
90% confidence interval. 

The lack of evidence for L x -T x evolution in 
our observation of MS 1054-0321 is consistent with 
the other evidence for Qo < 1* but test for 
Q 0 is not as powerful as the temperature function 


test, as discussed in Section 3.2. Some (luster 
evolution models predict little evolution in the Lx- 
Tx relationship in an Qo = 1 universe, as long as 
the minimum entropy of the X-rav gas remains the 
same for all clusters regardless of redshift (Kaiser 
1991, Evrard & Henry 1991, Bower 1997). The lack 
of evolution in Lx~T x would then require some sort 
of non-gravitational heating of the X-ray gas. 

4. SUMMARY 

The ASCA observations we have presented here 
imply that the EMSS cluster MS1054-0321 at 
z = 0.83 has an X-ray temperature of 12.3to 1 keV. 
Supporting optical redshifts of 12 member galaxies 
(this paper) and the weak lensing map (Luppino 
& Kaiser 1997) provide a velocity dispersion and a 
lensing mass that agree with the mass inferred from 
the high X-ray temperature. Our ROSAT image 
exhibits significant substructure, but numerical 
modelling of cluster evolution has shown that the 
presence of substructure does not preclude using a 
cluster’s X-ray temperature to estimate its virial 
mass (e.g., Evrard et al. 1996). 

This high a temperature in so distant a cluster 
is powerful evidence that Qo < 1* Using the 
integral Press-Schechter relation and the current 
temperature function of clusters, we calculate the 
expected number density of such hot clusters at 
z = 0.83, under the assumptions of Gaussian 
perturbations and Qo = 1- We find that the product 
of the expected number density and the effective 
volume of the EMSS in the z = 0.5 — 0.9 redshift bin 
is less than a few times 10” 5 . Similar calculations 
for two other hot EMSS clusters at z ^ 0.54 give 
probabilities less than a few times 10” 3 for each one. 

For Qo ~ 0.3, the temperature function of hot 
clusters is expected to evolve very little to redshifts 
> 0.5 (Eke et al. 1996), so that the hot EMSS clusters 
are quite consistent with a low-Q 0 universe. We also 
find that MS1054-0321 luminosity and temperature 
are consistent with the low-z Lx~Tx relation, a 
result also consistent with Qo < 1. These results 
strongly support the mounting body of evidence 
(e.g., Carlberg et al. 1997) showing that the mean 
density of matter in the universe is Qo ~ 0.3. 
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TABLE 1. ASCA OBSERVING INFORMATION 


Detector 

Aperture Radius 
Arcmin 

Count Rates 
s _1 

Useful Exposure Time 
s 


2.53 


62,019 


3.25 


58,176 


5.00 

0.012 

65,534 

GIS3 

5.00 

0.016 

65,920 


TABLE 2. MEMBER GALAXY REDSHIFTS 


ID 

Redshift 

Offset E/W (+/-”) 

Offset N/S (+/-”) 

l(cD) 

0.8309 ± 0.0008 

0.0 

0.0 

2 

0.8342 ± 0.003 

+23 

-6 

3 

0.8127 ± 0.0003 

+56 

-5.5 

4 

0.8213 ± 0.0007 

+37 

+30 

5 

0.830 : ±0.008 

-14 

-20 

6 

0.8209 ± 0.001 

-29 

-13 

7 

0.8286 ± 0.001 

-31.5 

-11.5 

8 

0.8353 ± 0.0006 

-38 

-8 

9 

0.8332 ± 0.001 

-43 

-6 

10 

0.8143 ± 0.002 

-38 

+0.5 

11 

0.8378 ± 0.003 

-80 

-44 

12 

0.8319 ± 0.002 

-97.5 

-38.5 
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disbursements less program income NET DISBURSEMENTS 

received) FOR REPORTING 


5. TOTALS (Should correspond with amounts shown on SF 272 as follows: column (c) the 
same as line 1 1 h; column (d) the sum of lines 1 lh and Hi of this SF 272 and cumulative 
disbursements shown on last report. Attach explanation of any differences.) 


272 -: 


PnMcnbod by Offic* of Min*t«m«vt and Budget 
Cir No. A-110 


EQUIPMENT AND PROPERTY INVENTORY REPORT 




*The cost of these was split 39%/61 % between this grant and NAG5-2615 



